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ABSTRACT 

The first Herschel Hi- Gal images of the galactic plane unveil the far-infrared diffuse emission of 
the interstellar medium with an unprecedented angular resolution and sensitivity. In this paper, we 
present the first analysis of these data in combination with that of Spitzer GLIMPSE & MIPS GAL. 
We selected a relatively diffuse and low excitation region of the I ~ 59° Hi- Gal Science Demonstration 
Phase field to perform a pixel by pixel fitting of the 8 to 500 jam SED using the DustEM dust emission 
model. We derived maps of the Very Small Grains (VSG) and Polycyclic Aromatic Hydrocarbons 
(PAH) abundances from the model. Our analysis allows us to illustrate that the Aromatic Infrared 
Bands (AIB) intensity does not trace necessarily the PAH abundance but rather the product of 
"abundance x column density x intensity of the exciting radiation field" . We show that the spatial 
structure of PACS 70 jam map resembles the shorter wavelengths (e. g. IRAC 8 jam) maps, because 
they trace both the intensity of exciting radiation field and column density. We also show that the 
modeled VSG contribution to PACS 70 jam (PACS 160 jam) band intensity can be up to 50% (7%). 
The interpretation of diffuse emission spectra at these wavelengths must take stochastically heated 
particles into account. Finally, this preliminary study emphasizes the potential of analyzing the full 
dust SED sampled by Herschel and Spitzer data, with a physical dust model (DustEM) to reach the 
properties of the dust at simultaneously large and small scales. 
Subject headings: methods: data analysis - ISM: dust - infrared: ISM 



1. INTRODUCTION 



The first Herschel (Pil bratt et all |2010) images of 
the galactic plane obtained during the Science Demon- 
stration Phase (SDP) reveal with an unprecedented 
beauty and detail the far-infrared diffuse emission of 
the interstellar med ium (ISM). The Hi-Gal project 
(Mol inari et al.ll201Q[ ) provides an unbiased photometric 
survey of the inner Galactic plane emission between 70 
and 500 jam where the spectral energy distribution (SED) 
is dominated by the largest dust grai ns (a ~ 0.1 jam). 
Comb ined with the Spit zer GLIMPSE (IChurchwell et"atl 
[2001 and MIPSGAL (jCarev et all I2009D surveys that 
cover the dust SED from 3.6 to 24 /im, light mostly emit- 
ted by the smallest grains (a< 10 nm), we will for the 
first time reliably sample the full dust SED of the diffuse 
emission at all spatial scales down to ~ 40" and over a 
large fraction of the galactic plane. 

At 70 and 160 jam many of the observed structures 

* HERSCHEL IS AN ESA SPACE OBSERVATORY WITH 
SCIENCE INSTRUMENTS PROVIDED BY EUROPEAN-LED 
PRINCIPAL INVESTIGATOR CONSORTIA AND WITH IM- 
PORTANT PARTICIPATION FROM NASA. SPITZER SPACE 
TELESCOPE IS OPERATED BY THE JET PROPULSION 
LABORATORY, CALIFORNIA INSTITUTE OF TECHNOL- 
OGY UNDER A CONTRACT WITH NASA 
Electronic address: compiegne@cita.utoronto.ca 

1 Canadian Institute for Theoretical Astrophysics, University of 
Toronto, 60 St. George Street, Toronto, ON M5S 3H8, Canada 

2 Spitzer Science Center, California Institute of Technology, 1200 
East California Boulevard, MC 220-6, Pasadena, CA 91125, USA 

3 Department of Astronomy Sz Astrophysics, University of 
Toronto, 50 St. George Street, Toronto, ON M5S 3H4, Canada 

4 Centre d'Etude Spatiale des Rayonnements, CNRS et Univer- 
site Paul Sabatier-Toulouse 3, Observatoire Midi-Pyrenees, 9 Av. 
du Colonel Roche, 31028 Toulouse Cedex 04, France 

5 INAF-Istituto Fisica Spazio Interplanetario, Via Fosso del Cav- 
aliere 100, 1-00133 Roma, Italy 



show a close spatial correlation with some of the classical 
infrared tracers of small dust particles, like the emission 
measured at 8 and 24 /im, and therefore, one expects 
a partially similar origin and/or a similar dependence to 
ISM physical properties. To properly understand the ori- 
gin of this emission is the first step of the data analysis. 

Dust is heated by stellar ultraviolet (UV)- visible pho- 
tons and re-radiate the absorbed energy in the infrared. 
Assuming some dust properties, this dust emission is 
widely used to trace crucial quantities like cloud masses 
or star forming activity. However, to derive accurately 
these quantities requires to take into account the dust 
properties evolution, how it affects its emission and how 
it reflects the ISM properties. On the other hand, dust 
plays a crucial role for the ISM physics and chemistry. 
Since dust and the rest of ISM are tightly interlock, 
changes on the dust properties affect those of the gas 
phase. It is then necessary to characterize the physi- 
cal processes responsible for the dust evolution to assess 
their role within the ISM lifecycle and to be able to use 
dust as a reliable tracer of the ISM. Accessing the full 
dust SED over a large fraction of the sky allows to follow 
the behavior of every dust component over a broad range 
of physical conditions which is an invaluable information 
for the study of the dust evolution. 

In this paper, we present the first analysis of the Hi- 
Gal data combined with that of GLIMPSE and MIPS- 
GAL to obtain the full dust SED from 8 to 500 jam. In 
the section[2j we describe the observations and the fitting 
method that makes use of a physical dust model, DustEM 6 
(jCompiegne et al.ll2010f ) . In the section[3j we present and 
discuss the results and how our analysis method allows us 

^ To be found at http://www.ias.u-psud.fr/DUSTEM 
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Fig. 1. — IRAC8/im, PACS 70 jum and SPIRE 500/im maps. The obtained maps for the fitted parameters are Ypah, YvSG, Nh and 
Ummp83- Also shown is a map of the product Ypah x Nh x ^MMP83-The two boxes delineate the areas used to obtain the spectra of 
Fig. [2] Note that the black spots are due to pixel masking (e.g. 8 /im point sources). 



to better interpret the behavior of extended emission ob- 
served with the Herschel/Spitzer photometers. We con- 
clude in the section|4j 

2. OBSERVATIONS & ANALYSIS METHOD 

Two fields were observed (at I ~ 30° and I ~ 59°) dur- 
ing the SDP as part of the Hi-Gal progra m with PACS 70 
and 160 jam and all the SPIRE channels (Moli nari et al.l 
l2010f ). For this study, to satisfy the assumptions made 
in the model we describe below, we focus on a subfield 
of the I ~ 59° field that shows no obvious HII regions 
or young stars cluster and that is relatively diffuse re- 
garding the rest of the I ~ 59° field and the I ~ 30° 
field. We also avoid the higher latitudes (|b| > 0.8°) that 
display weak intensity where Herschel data could be less 
reliable at this early stage of data processing. The maps 
were obtained with the ROMAGAL pipeline (Traficante 
et al., 2010, MNRAS, submitted). For PACS 160 /im, 
SPIRE 250, 350 and 500/im, we applied the calibration 
described in lBernard et al.l (j201Q[ ). The gain uncertainty 
for these data is taken to be 20%. The zero level was cor- 
rected by cross-calibration with Planck data with abso- 
lute uncertainties of ± 19.8, 14.6, 7.5 and 3.0 MJy sr _1 at 
160, 250, 350 and 500 /im, respectively, which represent 
about half of the gain uncertainty for the faintest pixels 
of the studied region. We cross-calibrate the PACS 70 /im 
data with the MIPS 70 /im data from the MIPSGAL sur- 
vey (Paladini et al., 2010, in prep). We then apply 
the 15% gain uncertainty of MIPS 70 jam. We also use 



the IRAC 8 jam data from the GLIMPSE survey and the 
MIPS 24 /am data from the MIPSGAL survey. The gain 
uncertainty for these two dataset is taken to be 10%. 
Point sources are subtracted from the IRAC 8 /im data. 

The zodiacal light c ontribution, import ant at shorter 
infrared wavelengths (jKelsall et al.l Il998h has been re- 
moved from all the Spitzer data (on average, 1.1, 16.5 
and 4.5MJysr _1 at 8, 24 and 70 /im, respectively) and 
then from the PACS 70 /im following its cross-calibration 
on MIPS 70 jam. We do not perform any subtraction of 
the zodiacal emission at longer wavelenghts where its 
contribution is negligible (<lMJysr _1 ). We bring ev- 
ery map to the lowest resolution of the SPIRE 500/im 
one (FWHM ^ 37 // ). We do this assuming a Gaus- 
sian Point Spread Function of appropriate width. We 
then project all maps into the SPIRE 500 jam grid (pixel 
field of view- 11.5"). Fig. □ displays the IRAC 8 /im, 
PACS 70 jam and SPIRE 500 jam maps. 

We use the DustEM dust model described in 
iCompiegne et al] (j2010h to analyse the data. We con- 
solidated the standard four grain populations into three: 
the Polycyclic Aromatic Hydrocarbons (PAH), Small 
Amorphous Carbons representing the Very Small Grains 
(VSG) and we merged Large Amorphous Carbons and 
Silicates in to a single Big Gr ains (BG) population. Using 
the mpfit ([Markwardt 2009) IDL minimization routine 7 , 
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Fig. 2. — Mean spectra over the boxes shown on Fig.[T]and cor- 
responding fitted model. The upper (lower) spectrum correspond 
to the northern (southern) box. For comparison, the blue lines 
show the spectra obtain for the fitted Njj and Ummp83 but for 
the reference diffuse high galactic latitude dust properties (DHGL, 
Ypah and Y VSG = 1). 

we choose to adjust the following four parameters to fit 
the observed SED for each pixel: the (i) PAH and (ii) 
VSG abundances relative to BG, Ypah and Yvsg, (hi) 
the BG opacity, tbg, and (i y ) Umm P83 a s caling factor 
of the "solar neighborhood" Mat his et al.l (jl983[ ) (here- 
after MMP83) exciting radiation field. Self-extinction 
along the line of sight can be significant at 8 jiva and is 

accounted for assuming I\ = Jo, a 1 ~^ a a where Jo, a is 
the integrated emissivity and t\ the total dust opacity 
that is computed from the dust model using Ypah, Yvsg 
and tbg- Since in this first analysis we do not perform 
any separation of the components along the line of sight, 
the derived parameters result from the spatial mixing of 
different physical conditions. Fortunately, at I ~ 59°, we 
expect only contributions from the Vulpecula star forma- 
tion region (d~2kpc) and the Perseus arm (d~8.5kpc). 

We focus on the smallest particles behavior (PAH and 
VSG) and hence, we do not assume any variation of the 
BG properties (i.e. tbg/Nh = constant). In that case 
and to ease the following discussion, we can convert tbg 
into the hydrogen column density, Nh, and assume the 
relative abundances of PAH and VSG, Ypah and Yvsg 
to be their abundance relative to hydrogen. Indeed, emis- 
sivity and /or abundance of the BG is know n to evolve 
(e.g. iStepnik et all 120031 iDesert et a_Lll2008h and can be 
studied using the Herschel data (see lParadis et al.ll2010[ ) 
but to take such variations into account is beyond the 
scope of this paper. 

Fig. [2] illustrates the obtained fitted spectra for two ob- 
served SEDs. The photometric points computed from 
a modeled spectrum take into account the color correc- 
tions. Schematically, the 160, 250, 350 and 500 fim pho- 



tometric points give constraints on Ummp83 and tbg 
(Nh) through the BG emission. The shape of the BG 
SED depends on Ummp83 as this grain population is 
at thermal equilibrium, in contrast with the stochasti- 
cally heated grains (PAH and VSG) whose SED shape 
is invariant regarding the radiation field intensity (see 
iDraine fc Lill2007l : ICompiegne et al.ll201Q[ ). Hence, given 
Ummp83 the absolute level of the BG SED constrains 
Nh- Finally, the 8 and 24 /am constrain the abundances 
of the two stochastically heated population since the in- 
tensity of their emission scales as Y x Nh x Ummp83- 

3. RESULTS & DISCUSSION 

Fig. Q] shows maps of the Ypah, Yvsg, Nh and 
Ummp83 parameters and a map of the product Ypah x 

N H X U M MP83- Xleduced ~ 2 for a11 P ixels over tnese 

maps. Nh is given in unit of 10 20 H cm~ 2 and Ummp83 
is dimensionless. Ypah and Yvsg are given rela- 
tive to the value for the diffuse high galactic latitude 
medium (DHGL, \b\ > 15°), M PA h/M h = 7.8 10~ 4 an d 
M V sg/Mh = 1.65 10~ 4 fsee ICompiegne et all (201(1 ). 
Notice that M B g/M h = 9.25 10" 3 . The uncertainties 
on these parameters are given as computed in mpfit from 
the covariance matrix. 

The column density lies between ~1.5 10 22 H cm~ 2 for 
the most diffuse part and ~1.1 10 23 H cm~ 2 toward the 
dense filamentary structures, with a relative uncertainty 
of ~18%. The Ummp83 ranges between ~0.5 and ~1.5 
with a relative uncertainty of ~22%, which is consis- 
tent with the absence of young stellar cluster over our 
field. In agreement with the fact that the radiation 
field is shielded, we see that Ummp83 decreases toward 
the highest column density regions. These Ummp83 and 
Nh values are in a greement with previous estimates (see 
iBernard et al.ll2010l ). Ypah goes from ~ 1.0 to ~2.0 with 
a relative uncertainty of ~16%. The obtained Yvsg 
spans a wider range than Ypah-, going from ~0.5 to ~3.0 
but with a relative uncertainty going from ~40% for the 
highest values to ~60% for the smallest values. The vari- 
ations of these parameters over the map are therefore 
significant. 

Fig. [2] shows the averaged SED over the two boxes 
seen on Fig.[T] and illustrate the wide Yvsg variations. 
For the top (bottom) panel spectrum we have Nh = 
5.1±0.9 10 22 # cm- 2 (2.5±0.5 10 22 i^cm- 2 ), U M MP83 = 
1.2±0.2 (0.9±0.2), Ypah = 1.7±0.2 (1.4±0.2) and 
Yvsg = 0.8±0.5 (2.7±1.0). 

Ypah jYysG seems to vary at large spatial scales (de- 
creases from b ~ 0.2° to b ~ —0.7°) and also at the 
edge of some dense filaments (near the "chimney" region 
around coordinate 58.4, +0.55). Both Ypah and Yvsg 
decrease towards some of the densest filamentary struc- 
tures but it appears not to be systematic and could be 
biased by the too simplistic assumption made on the self- 
absorption. A decrease of the smallest dust abundance 
towards dense regions regarding the biggest grains was 
reported by previous work and was interpreted as the 
coagulation of th ese smallest particles together with th e 
bigger ones fe.g. lStepnik~eTaI . 2003; Flage v et alJ l2009). 
Previous studies also reported a lack of correlation or 
even an anti-correlation between the aromatic infrared 
bands and the mid-IR continuum emission (interpreted 
as the evolution of small dust properties) at the illumi- 
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nated ridge of mo l ecular clouds (e.g. Abergel et al.ll2002t 
iBerne et all 120071 : iCompiegne et al .1 l2008h . toward high 
galactic cirrus at the interface between atomi c and molec- 
ular material (jMiville-Deschenes et al.ll2QQ2l) or at galac- 
tic sc ale in the Large Magellanic Cloud (|Paradis et al.l 
l2009h . 

The spatial structure of SPIRE 500 jam intensity, that 
is dominated by the BG contribution, is well correlated 
with Nh- Indeed, Ummp83 and subsequently the BG 
temperature is quite stable in our field so that the in- 
tensity variations at these wavelengths (Rayleigh tail of 
the blackbody like emission) are dominated by the col- 
umn density variations (assuming no BG emissivity vari- 
ations). Comparing the 8 /am and Ypah maps, it is strik- 
ing that the AIB intensity does not directly trace the 
PAH abundance. PAH are stochastically heated and the 
observed AIB intensity (i.e. IRAC 8 jam) then scales lin- 
early with the product Ypah x Nh x Ummp83 as seen 
on Fig.[TJ The only difference is related to the extinction 
along the line of sight that can be seen in the IRAC 8 jam 
(as dark filamentary structures correlated with the Nh 
map) and not in the Ypah x Nh x Ummp83 map. The 
PACS 70 jam image shows a close spatial correlation with 
the classical infrared tracers of small dust particles, like 
the emission measured at 8 /im (see Fig. [I]). The emis- 
sion at 70 /im is not due to a single dust component. 
As seen on Fig[2j in the framework of our model both 
the VSG that are stochastically heated and the BG that 
are at thermal equilibrium contribute to the emission at 
this wavelength. For Ummp83 ^ 100, the BG emission at 
70 jam fall in the Wien part of the blackbody like emission 
that makes it more sensitive to Ummp83 than at 500 jam. 
Therefore, the two emission components (from VSG and 
BG) at 70 /im are sensitive to both Nh and Ummp83 ex- 
plaining the better correlation of this map with the 8 /am 
map (also sensitive to both Nh and Ummp83) than with 
the 500 jam map (more sensitive to Nh)- 

The model allows us to compute the relative contribu- 
tion of the three dust populations in the different pho- 
tometric bands. The VSG contribution to PACS 70 jam 
increases if the VSG abundance relative to BG increases 
and/or if Ummp83 decreases (shifting the BG emission 
toward longer wavelengths). In the studied field, the 
VSG contribution to the PACS 70 jam intensity goes from 
~ 10% up to ~ 50% with a median value of ~ 27%. 
The maximum contribution of VSG to PACS 100 jam and 
PACS 160 jam is - 17% and - 7% (the median is 9% 
and 3%), respectively. For the top (bottom) panel spec- 
trum of Fig.EJ the contribution is 12, 4, 2% (35, 13, 
5%) for PACS 70, 100 and 160 jam, respectively. This 



result strongly suggests that the proper analysis of Her- 
schel spectrum of diffuse emission including 70 jam may 
require to account for stochastically heated grains. 



4. SUMMARY & CONCLUSION 

We have presented the first analysis of the diffuse emis- 
sion of the galactic plane as observed by Herschel com- 
bining the Hi-Gal data with the GLIMPSE/MIPS GAL 
Spitzer data. Toward a subfield of the 1~ 59° Hi-Gal 
SDP field, we performed a pixel by pixel fitting of the 
full dust SED between 8 an d 500 /am using a physi cal 
dust emission model, DustEM (jCompiegne et al1l2010h . 

Assuming that the BG properties remain constant, the 
unique wavelength coverage provided by the Spitzer and 
Herschel photometric observations allows us to derive 
the following parameters for our dust model : the PAH 
and VSG abundances, Ypah and Yvsg, the column den- 
sity, Nh, and the intensity of the exciting radiation field, 
Ummp83- To our knowledge, this is the first time PAH 
and VSG abundance maps are derived using such a pixel 
by pixel SED fitting at resolution < V and over such an 
extended field. These abundances, as well as Nh and 
Ummp83, vary significantly over the field. As already re- 
ported by previous studies, Ypah and YysG appears not 
to be positively correlated. 

Although it was already theoretically known, our 
analysis method provide a firm demonstration that 
IRAC 8 /am does not trace the PAH abundance but the 
product Ypah x Nh x Ummp83- We also showed that 
at 70 /im, the modeled emission is due to both the BG 
and the VSG. At these wavelengths the BG emission is 
sensitive to both the intensity of the exciting radiation 
field and the column density likewise the VSG, explain- 
ing the similar spatial structure seen in PACS 70/im maps 
regarding shorter wavelengths. Using our model, we de- 
rived the VSG contribution to the PACS channels that 
can be up to - 50%, - 17% and - 7% at 70, 100 and 
160 /im, respectively. We conclude that the interpreta- 
tion of Herschel spectrum of the diffuse emission down to 
70 jam may require to take into account the stochastically 
heated population (VSG). 

Finally, our analysis allows for a better understand- 
ing of the first Herschel images of the galactic plan dif- 
fuse emission by disentangling between the different dust 
population contributions and also revealing the great po- 
tential of the Herschel/ Spitzer synergy combined with a 
phys ical dust model for the study of dust evolution (see 
also lAbergel et all l2010h . 
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